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Abstract 



Motivated by the recent experiment of Hochgrafe et al, we have investigated 
the magnetoplasmon excitations in a periodic array of quantum wires with a 
periodic modulation along the wire direction. The equilibrium and dynamic 
properties of the system are treated self-consistently within the Thomas- 
Fermi-Dirac-von Weizsacker approximation. A calculation of the dynamical 
response of the system to a far-infrared radiation field reveals a resonant 
anticrossing between the Kohn mode and a finite-wavevector longitudinal ex- 
citation which is induced by the density modulation along the wires. Our 
theoretical calculations are found to be in excellent agreement with experi- 
ment. 
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I. INTRODUCTION 



Modern microstructuring techniques have made it possible to experimentally study a 
wide variety of low- dimensional electron systems such as quantum dots, wires and rings.0 
These systems are usually fabricated starting from a two-dimensional electron gas (2DEG) 
located at a modulation-doped semiconductor heterojunction. The application of laterally 
microstructured metallic gates@ or the use of deep-mesa etching techniques! allow for the 
precise patterning of the 2DEG. Over the past 15 years, a great deal of information has been 
obtained about the transport and far-infrared (FIR) optical properties of these systems. For 
example, magneto-transport measurements of laterally modulated 2DEG's reveal interesting 
magnetoresistance oscillations associated with the commensurability of the cyclotron radius 
and modulation period a.0 Far-infrared transmission experiments, on the other hand, re- 
veal a rich spectrum of collective excitations whose frequencies exhibit novel magnetic field 
dispersions.a 

In this paper, we are interested in the magnetoplasmon excitations of an array of parallel 
quantum wires which are also periodically modulated along their length. Such a system was 
recently fabricated by Hochgrafe et a/.,i and studied using FIR transmission spectroscopy. 
They observe a pair of collective excitations which appear to interact and anticross as a 
function of magnetic field. This behaviour is qualitatively explained by elementary con- 
siderations of the magnetoplasmon mode spectrum of a uniform wire as described by the 
theory of Eliasson et a/H The periodic modulation along the length of the wires acts as a 
grating-coupler which couples the Kohn mode to a finite-wavevector ID magnetoplasmon of 
the uniform wire. The experimental data also show a novel kink-like feature in the vicinity 
of the cyclotron frequency uj c , that was neither discussed, nor explained in Ref. [5]. 

Our purpose here is to theoretically investigate the system studied by Hochgrafe et al. 
in order to better understand the nature of the magnetoplasmon excitations observed. In 
particular, we would like to gain some insight into the nature of the periodic modulation 
of the electronic density in these systems, and how this leads to the coupling of the modes 
discussed above. To this end, we adopt a hydrodynamic theory based on the Thomas- 
Fermi-Dirac-von Weizsacker (TFDW) approximation. This approach has already been used 
to describe the collective excitations in laterally modulated 2DEG'si as well as other low- 
dimensional electronic structures.iH0 Its main virtue is that it treats the ground state and 
dynamic properties of the system in a consistent fashion.ll! Our earlier successful applications 
give us confidence that the approach is able to accurately describe the collective excitations 
in the systems studied in Ref. [5]. 

The plan of the paper is as follows. In Sec. |H| we define our model for a periodically 
modulated quantum wire array. Then, we systematically analyze the effects of this modu- 
lation on the ground state, and dynamical properties of the system. To the extent possible, 
we make direct comparisons with the experiment of Ref. [5] . Finally, in Sec. |TT|, we present 
our concluding remarks. 

II. MODULATED WIRES 

The TFDW approximation to theground state of a quantum confined 2DEG is given by 
(atomic units are used throughout )0E~ 
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The first term is the Thomas-Fermi approximation to the kinetic energy, the second term 
is a von Weizsacker-like gradient correction, and the third term is the Dirac local exchange 
energy. For simplicity we neglect a correlation energy contribution which in any case has a 
small effect. The last two terms characterize the Hartree self-energy of the electrons and the 
interaction with the external potential, respectively. Once the confining potential v ext (r) has 
been defined, the variational minimum of Eq. (|ip yields the ground state electronic density, 
no(r), for the system of interest .100 To describe the array of modulated quantum wires in 
[5] we adopt a conning potential having the formO 

Vafi(x, y) = \ulx 2 + V y cos(Gy) , (2) 

with G = 2ix/a. It is understood that the potential in Eq. (||]) is centered in each unit cell and 
periodically extended throughout the array with period a. Setting V y = describes a periodic 
array of parabolically confined quantum wires, whereas V y ^ provides a modulation along 
the wire direction, also with period a (see Fig. |I|). Keeping in mind that the samples in 
[5] are prepared from a modulation-doped AlGaAs/GaAs heterostructure, the appropriate 
physical parameters are: the effective Bohr radius = 103 A and the effective Rydberg of 
energy, Ry* = e 2 /2eao = 5.4 meV. In these units, a 2D electronic density of 1.0 x 10 n cm -2 
is given by n 2 n = 0.10(o,q) _2 , and an energy of 0.5 a.u = lRy* is equivalent to 43.57 cm -1 . 

In Fig. |I| we present the equilibrium density distributions for a small section of the 
array. Figure |l|(a) is generated with the following physical parameters: an average density 
n2D = 1-0 x 10 11 cm -2 , a = 600 nm, cu = 0.427 a.u, and V y = 1 a.u. As will be discussed 
in more detail later, this choice allows us to reproduce much of the experimental resultsi 
obtained with the sample having a density of n 2 v = 1.2 x 10 11 cm -2 = 0.12(cIq) _2 . We 
therefore believe that the very strongly modulated density shown is typical of the densities 
studied in the experiments. In Fig. |T](b), we show the effect of keeping the confining potential 
fixed and doubling the 2D density to 7I2D = 0.2. It can be seen that the degree of corrugation 
is reduced in both directions, and that the densities from adjacent wires are now beginning to 
overlap. In going from Fig. |l|(b) to Fig. |l|(c) we have doubled the strength of the modulation 
to V y = 2 a.u. This gives rise to a much stronger modulation of the density along the length 
of the wires which now appear as a string of quantum dots. This is the density profile used to 
simulate the FIR response of the system as described in the second half of the experimental 
study in Ref. [5]. Finally, in Fig. [l|(d), we show the density profile generated with the same 
potential parameters as in Fig. |l|(c), but with a further 50% increase in the average density. 
This density takes on the character of an array of coupled quantum dots. 

The dynamics of a mesoscopic system within the TFDW approach is governed by a set 
of (linearized) hydrodynamic equations,000 

^+V'M = 0, (3) 
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with 8n(r, t) = n(r,t) — n (r), and 

<9v 

_ = -TV - VU mt + F ext . (4) 
dt w 

Eq. (|) is the usual continuity equation, and Eq. (f|) is a statement of Newton's second law. 
The quantity U mt plays the role of an effective potential for the internal force fluctuations, 
and F ext is the force due to any additional external electromagnetic fields imposed on the 
system, including a homogeneous, static magnetic field oriented normal to the 2DEG. In 
Eq. (f|), we have also included a phenomenological viscous damping term, —7V, to account 
for Drude scattering. This term has the effect of giving the collective modes a finite lifetime. 

The solution of Eqs. (HD and (W) allows us to calculate the physically relevant, time- 
averaged FIR power absorption.l'0Oll In terms of the induced current density j md (r,t) = 
— n (r)v(r, t), we have 

(P) t = Q d 2 r j ind (r,t) -E cxt (r,i)^ , (5) 

where 

E cxt (r, t) = (E x k + E y y) cos(arf) (6) 

is a spatially uniform radiation field polarized along (E x = 0) or perpendicular to (E y = 0) 
the direction of the wires. Most of our calculations are carried out for the perpendicular 
polarization, but some calculations are also done for the parallel polarization in order to 
identify some of the longitudinal excitations at low magnetic fields (B ^ 0.5T). 

The theoretical FIR power absorption for the density profile corresponding to Fig. [l|(a) 
is shown in Fig. The curves are displayed in steps of AB = 0.25T starting at B = 0T, 
and are offset vertically for clarity. The dashed lines are guides to the eye and indicate 
the peak trajectories; the open symbols at the end of each line serve to label the differ- 
ent magnetoplasmon modes (see Fig. |3|). The B = 0T trace shows an isolated peak at 
uj ~ 37 cm -1 which is very close to the parabolic confinement frequency uj . This is a clear 
signature of the center-of-mass Kohn mode (KM) predicted by the generalized Kohn theo- 
rem,@ and corresponds to a rigid oscillation of the equilibrium density perpendicular to the 
wires. With increasing magnetic field, this peak at first shifts slightly to higher frequencies, 
but around B ~ 0.8T, this trend is reversed with the appearance of a small shoulder to the 
high frequency side of the peak. The shoulder evolves into a separate peak which gains in 
oscillator strength and shifts to higher frequencies along the dashed line ending with the 
unfilled diamond. Although difficult to resolve, there is also a weak peak labeled by the 
unfilled triangle. 

The positions of the absorption peaks as a function of magnetic field are displayed in 
Fig. |3|. The unfilled symbols used in this figure correspond to the labels used to track the 
peak trajectories in Fig. |2]. For comparison, the filled squares are the raw experimental data 
taken from Ref. [5]. The dashed lines in Fig. |3| are the cyclotron frequency uj c and 2uj c . We 
will discuss the solid curves below. The splitting of the upper experimental mode between 
B = 1.5T and 2.5T is attributed in [5] to the presence of the Bernstein mode at uj = 2cu c ^^ 
Apart from this detail, it is clear that the TFDW results provide a smooth interpolation 
through the experimental points, and overall are in excellent agreement with the data. In 
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particular, the size and position of the anticrossing between the two lowest-lying branches 
(open circles and open diamonds) is well reproduced. 

In order to explain the details of the magnetic dispersions in Fig. [| it will prove useful 
to first consider the simpler case of a 2DEG which is modulated in one direction (see Ref. [7] 
for details). These systems were generated by setting ojq = in Eq. (|2|) and taking the 
modulation to be in the x-direction, i.e., V y — > V x . An additional weak density modulation 
along the wires will induce a coupling between the q y = modes of the wire to modes with 
wavevector q y = G = 2n/a. It is therefore of interest to see how the q y = 2n/a mode 
frequencies evolve as a function of magnetic field. In Figs, ^(a)-(c), we present the magnetic 
dispersions of a uniaxially modulated 2DEG at fixed q y = 2n/a. The excitation spectrum for 
the weakly modulated system (V x = 0.25 a. u) is shown in Fig. £|(a) and essentially exhibits 
the magnetoplasmon dispersion of a homogeneous 2DEG. The dense spectrum of modes 
starting from u = at B = 0T are a result of the broken translational invariance in the 
modulation direction, and eventually evolve into a series of edge magnetoplasmons (EMP's). 
In Fig. f|(b) we show the magneto-dispersion for the moderately modulated 2DEG having 
a modulation V x = 1.25 a. u. For this strength of modulation, one of the magnetoplasmon 
branches starting from u ~ 36 cm -1 acquires a weak negative dispersion and exhibits an 
anticrossing-like behaviour with a mode which has the character of an EMP.0 The mode 
dispersing upward from uj ~ 45 cm" 1 is associated with the KM and continues to increase in 
frequency with increasing magnetic fieldB Finally, in Fig. f|(c) we show the spectrum of the 
q y = 2ir/a modes for a modulation amplitude of V x = 1.5 a. u. In this modulation regime, 
the 1DMP has a pronounced negative dispersion for B < 2 T, and the gap associated with 
the anticrossing between the q y = 2n/a 1DMP and EMP modes has decreased. 

We are now in a better position to comment on the details of the magnetic dispersions 
in Fig. ^. Several solid curves are shown which correspond to the smooth wire limit. The 
solid curve which is lowest at B = is the dispersion relation of the q y = center-of-mass 
KM, namelyi'0 

uj 2 (B)=uj 2 +ujI (7) 

with ujq = 37.2 cm" 1 . The next two solid curves are the magnetic dispersions of the two 
q y = 2n/a modes calculated by taking V y — and u>o = 37.2 cm -1 in Eq. (g) and are 
analogous to the pair of modes starting from uj ~ 45 cm -1 in Fig. 4(c). We thus conclude 
that the two branches denoted by open triangles and diamonds in Fig. [5] are derived from 
the q y = 2tc /a KM and the q y = 2n ja 1DMP of the uniform wire array as discussed above. 
However, in the present situation, the q y = 2ir/a 1DMP is in fact coupled to the q y — KM 
by the periodic modulation along the wires, which results in the anticrossing behaviour seen 
in the vicinity of B = 1 T. For B < IT, the mode indicated by open circles has the character 
of the KM, while for B > IT, the mode changes over into the q y = 2n/a ID EMP. A similar 
change in character occurs for the open diamond branch, but in the opposite sense. This is 
precisely the interpretation of the experimental results given in Ref. [5] . The highest energy 
branch (open triangles) corresponding to the weak high energy resonance in Fig. |3] is clearly 
the q y = 2n/a mode belonging to the upward dispersing bulk mode (see e.g., Fig. |](c)). 
Finally, we note that the track of the lowest mode in Fig. |3] shows a small 'kink' at uj = u c . 
As we shall see below, this kink becomes more prominent as the density of the system is 
increased. We defer further discussion of this feature until then. 
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We next consider the effect of increasing the electronic density. This is achieved exper- 
imentally by means of illumination, and in one particular run resulted in an approximate 
doubling of the average density. To simulate this situation, we considered a system having 
a density of = 0.2 and the same confining potential used to generate the results in 
Fig. [| The density distribution in this case is shown in Fig. |l](b) where it can be seen that 
the modulation of the density along the wires is weaker than for 7X2D = 0.1. However, the 
simple assumption of an increase in the electron density was not able to account for the 
observed trends in the experimental data. Indeed, agreement with experiment could only be 
achieved by simultaneously increasing the amplitude of the modulation potential to V y = 2 
a. u. The corresponding density distribution is shown in Fig. |l](c). The modulation of the 
density along the wires is clearly more pronounced than in Figs. |l](a) or |l](b). In Fig. ||, 
we present the calculated FIR power absorption. As in Fig. |2|, there is a strong KM peak 
at B = 0T, consistent with the parabolic confinement in the x-direction. In addition, due 
to the increased coupling of the density in a direction perpendicular to the wires, we pick 
up some weak Drude absorption at cu — 0. The dashed curve labeled by the open square 
symbol tracks the location of this peak which follows closely the dispersion of the cyclotron 
mode to = u c up to the point where it interacts weakly with the 1DMP (open circles). 

The peak positions as a function of B are shown in Fig. [(]. As compared to Fig. |3|, 
the anticrossing between the KM and the 1DMP has now been pushed down in magnetic 
field so that there is already a strong interaction at B = T. As a result, the frequency 
of the lower mode indicated by open circles decreases monotonically for the entire range of 
magnetic fields presented, and at higher magnetic fields has the character of the ID EMP. 
The mode indicated by open diamonds attains a noticeable signal at B pa 0.8T as it evolves 
into the KM. However, unlike the case of lower electron density, the oscillator strength of the 
KM is shared with higher lying modes indicated by the open triangles and stars in Fig. 
Although the triangle mode in the figure is very weak, it appears more prominently in the 
FIR absorption spectrum when the radiation is polarized along the wires. This mode is the 
same q y = 2n/a mode discussed earlier (triangles in Fig. |3|) but in this case it interacts with 
another mode near B = 1.5T and continues along the starred branch at higher magnetic 
fields where its strength is comparable to that of the KM. 

The agreement with the experimental data (filled squares) is again very good, except in 
the vicinity of the Bernstein mode at 2u c . Our model provides a smooth and continuous 
interpolation through this region. The lower experimental branch closely follows the 1DMP 
indicated by open circles while the upper branch is clearly associated with the KM (open 
diamonds) . 

The weak structure in the trajectory of the lower experimental branch seen near u c 
can be explained in terms of the pair of q y = 2ir/a modes which appear to anticross in 
Fig. 4(c). In the present situation, these modes are excited by virtue of the 2ir/a modulation 
along the wire length. To confirm our interpretation of this kink-like feature, we considered 
a further increase in the density of the system without changing the potential used to 
generate Fig. |I|(c). The density profile associated with this situation is shown in Fig. |T](d). 
Figure [7| illustrates the FIR power absorption of the system, and Fig. || shows the peak 
positions as a function of the magnetic field. Comparing Fig. || with Fig. f|(c), we see that 
these dispersions share several common characteristics. In particular, we see that the open 
diamond and triangles in Fig. |3| have the same qualitative behaviour as the two q y = 2n/a 
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modes starting at uo ~ 45 cm -1 in Fig. |4](c). Also, the strong anticrossing between the open 
circles and open squares in Fig. || clearly has the same origin as the anticrossing between 
the 1DMP and the EMP mode in Fig. |](c). Thus, the weak "kink" in the experimental 
data of Ref. [5] is associated with the anticrossing of the q y = 2ir/a 1DMP and EMP modes 
which are excited by the G = 2n/a modulation along the wire direction. As Fig. 4 indicates, 
the apparent coupling between the 1DMP and EMP modes decreases as the magnitude of 
the uniaxial modulation increases. This explains why the "kink" feature is more apparent 
at higher densities and gets progressively weaker as the wires become increasingly isolated 
from each other. 

III. CONCLUSIONS 

In this paper, we have studied the magnetoplasmon excitations in an array of periodically 
modulated quantum wires within the TFDW approximation. We have shown that the 
modulation along the wire direction effectively acts as a grating coupler which induces a 
coupling between the q = Kohn mode and the q = (0,27r/a) ID plasmon, resulting in 
a resonant anticrossing of these modes at a finite magnetic field. We have compared our 
theoretical results with the data of Ref. [5], and found excellent agreement for all density 
and magnetic field regimes investigated. In addition, we have provided a simple explanation 
for the kink feature at uj = u c , which was observed, but not explained, in Ref. [5]. Based 
on the overall agreement between theory and experiment, we believe that the underlying 
ground state densities used in the calculations, as illustrated in Fig. 1, are in fact a realistic 
representation of the density profiles produced in the experiments. 
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FIGURES 



FIG. 1. Equilibrium density distributions, presented as surface and contour plots, for: (a) 
^2D = 0.1, V y = 1 a.u, (b) n 2 rj = 0.2, V y = 1 a.u, (c) n 2 r> = 0.2, V y = 2 a.u, and (d) 1%2d = 0.3, 
V y = 2 a.u. In all cases, a = 600 nm and lu = 0.427 a.u. The parameters pertaining to v cxt were 
chosen to reproduce the B = 0T excitation frequencies in Ref. [5]. 

FIG. 2. Calculated FIR power absorption of the modulated quantum wire array corresponding 
to the density in Fig. [j](a). The unfilled symbols are used to label the various peak trajectories, 
denoted by the dashed lines. The incident radiation is polarized in a direction perpendicular to 
the wires. 

FIG. 3. Theoretical B dispersion of the resonant FIR absorption data in Fig. ^ The unfilled 
symbols correspond to those in Fig. ||, and the filled squares are the experimental data from Ref. [5] . 
The dashed and full curves are explained in the text. 

FIG. 4. Calculated B dispersions of the modes in a uniaxially modulated 2DEG at a fixed 
wavevector q = (0, 27r/a). The amplitudes of the modulating potentials are (a) V x = 0.25 a.u. , 
(b) V x = 1.25 a.u. , and (c) V x = 1.5 a.u. The dashed line represents the cyclotron frequency. 

FIG. 5. As in Fig. || but for an equilibrium density defined by Fig. |l|(c). 

FIG. 6. As in Fig. |3|, but with the data extracted from Fig. |5[ 

FIG. 7. As in Fig. [2|, but for an equilibrium density defined by Fig. 0(d). 

FIG. 8. As in Fig. but with data extracted from Fig. 0. 
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